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FOREWORD 

A  cylindrical  coordinates  approach  to  phenomenological  band  structure 
calculations  introduced  previously  is  generalized  to  include  spin  degrees  of 
freedom  and  the  spin  orbit  interaction.  The  pertaining  reduction  of  the 
Schrodinger  equation  to  a  set  of  1-d  wave  equations  is  derived.  For  the  test 
case  of  PbSe  band  along  the  F-L  line,  the  correct  bands  are  reproduced  in  the 
numerically- stable  range  of  parameters. 
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discussions.  This  work  was  supported  by  the  Office  of  Naval  Research  and  the 
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ABSTRACT 

A  rapidly  convergent  method  for  band  structure  calculations,  based  on  a 
cylindrical  coordinates  expansion,  is  generalized  to  include  the  spin-orbit 
interaction.  This  approach  is  advantageous  particularly  for  materials  highly 
anisotropic  in  one  direction.  The  pertinent  wave  equations  are  tested  for  the 
PbSe  low  bands  along  the  F-L  line.  The  results  compare  well  with  PbSe  known 
band  structure  except  where  the  bands  group  together  to  form  degenerate,  or 
near  degenerate,  equal  parity  clusters. 
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SECTION  1 
INTRODUCTION 


In  a  previous  paper^  (hereafter  referred  to  as  I),  we  introduced  a 
method  for  band  structure  calculations  along  a  high- symmetry  axis.  The  method 
constitutes  of  expanding  the  wave  function  (and  potential)  in  cylindrical 
coordinate  multipole  and  reducing  the  three-dimensional  Schrodinger  equation 
to  a  set  of  coupled  one-dimensional  wave  equations  for  these  multipoles.  The 
primary  advantage  of  the  method  is  its  good  convergence,  in  particular  for 
anisotropic  systems  with  a  unit  cell  largely  extending  in  one  direction,  z'. 
Examples  for  such  systems  are  superlattices ,  intercalates,  and  a  lattice 
bounded  by  a  surface.  The  method  applies  equally  to  bulk  band  calculations 
along  a  high-symmetry  axis.  In  this  paper  we  generalize  the  method  to  include 
spin  degrees  of  freedom.  Such  an  extension  is  necessary  in  situations  where 
the  spin-orbit  (SO)  interaction  is  important.  A  case  in  point,  the  heavy- 
rlement  narrow  gap  semiconductors  belonging  to  the  II-VI  and  IV-VI  systems.^ 
For  these  materials,  the  SO  interaction  is  essential  to  correctly  obtain  the 
band  jap.^  Another  application  is  the  calculation  of  g-f actors  where  the 
inclusion  of  spin  degrees  of  freedom  is  obviously  indispensable.  We  outline 
first  the  method  without  spin  and  point  out  the  modifications  mandated  by  the 
presence  of  the  spin  degrees  of  freedom.  We  then  comment  on  the  main  points 
borne  out  by  the  band  calculation  of  PbSe,  which  serves  as  a  test  case. 


At  the  core  of  the  method  is  the  cylindrical  multipole  expansion  of 
lattice  periodic  functions.^  Examples  of  such  functions  are  band  wave 
functions,  the  pseudopotential,  and  the  spin-orbit  form  factor.  These 
multipoles  follow  naturally  from  the  introduction  of  cylindrical  coordinates, 
chosen  such  that  its  z-axis  coincides  with  z',  the  high-symmetry  axis.  For 
example,  the  multipole  expansion  of  a  band  wave  function  #(r')  where  the 
crystal's  momentum  Ic  points  at  the  z' -direction  is^ 


’l'(r')  -  I  I  't^(gp,z')  J^(gpP)  e^'*'**  (1-1) 

gp 

In  Equation  (1-1),  gp  runs  over  a  series  of  lattice-specific  nonnegative 
values  (Table  1-1),  J^(x)  is  the  Bessel  function  of  order  f  and  (p,4>.z')  are 
the  standard  cylindrical  coordinates.  Rapid  convergence  of  Equation  (1-1) 
implies  that  a  small  number  of  distinct  raultipoles  ♦^(gp.z')  yields  a  good 
description  of  4'(r'). 
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TABLE  1-1.  STRUCTURE  OF  fee  g-HEXAGONS  ASSOCIATED  WITH  THE  [111]  DIRECTION 

g  -  f/a  (in  A'^) 


f 

g(PbSe)(^^ 

(degrees) 

1. 

0. 

0. 

0. 

2. 

10.2606  - 

2,J| 

1.6754 

0. 

3. 

17.7714  = 

2.9018 

30 

4. 

20.5207  - 

-IF 

3.3507 

0. 

5. 

27.1465 

4.4326 

19.1066 

6. 

30.7813 

5.0261 

0. 

7. 

35.5430 

5.8036 

30 

8. 

36.9944 

6.0406 

13.8979 

9. 

41.0414 

6.7014 

0. 

10. 

44.7239 

7.3027 

36.5868 

11. 

47.0193 

7.6775 

10.8934 

(a)  We  used^  a(PbSe,  T  -  300K)  -  6.1243  A 

NOTE;  The  g  values  and  <^Q(g)  angle  (Equation  (3-1))  for  the  PbSe  (fee)  lattiee 
in  the  [111]  direetion.^ 
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The  convergence  of  Equation  (1-1)  hinges  on  two  general  properties, 
which  severely  limits  the  number  of  contributing  f  and  gp  indices.  The  £ 
variation  is  controlled  by  the  symmetry  of  the  z'-axis.  For  example,  group- 
theoretic  considerations  imply  that  for  a  fee  lattice  and  z'||[lll]  only  two  i 
values  need  be  considered.^  All  other  multipoles  are  phase  related  to  those 
two  multipoles.  As  a  consequence,  expansion  (1-1)  is  in  effect  an  expansion 
in  gp.  The  gp  expansion,  in  turn,  is  controlled  by  an  energy-geometry 
argument.  To  see  this  point,  note  chat  in  Equation  (1-1)  the  transversal 
behavior  of  ^'(r')  is  controlled  by  the  oscillations  of  J^(gpp),  with  period 
gp^.^  Therefore  for  "low"  band  energies,  the  important  terms  in  Equation  (1- 
1)  involve  only  the  lowest  gp- values,  up  to  a  cutoff  value  g^.  Based  on  the 
behavior  of  the  J^(x)  curve,  the  latter  has  been  estimated  in  I  to  satisfy  (i 
+  3/2)/g(-  =:  .5  a-j.  where  i  is  a  typical  i -value  and  is  of  the  order  of  the 
lattice  constant  in  the  transversal  direction.  For  the  PbSe  example  (see  I), 
g^  turns  to  be  approximately  the  second  g  value  (see  Table  1-1;  the  first  g 
value,  g  =  0,  always  contributes).  Thus  both  the  i  and  g  variation  in 
expansion  Equation  (1-1)  are  severely  restricted,  hence  the  rapid  convergence. 
It  is  important  to  emphasize  that  the  above  considerations  are  valid 
regardless  of  the  extension  of  the  unit  cell  in  the  z'  direction!  In 
particular,  the  good  convergence  of  Equation  (1-1)  applies  equally  for  bulk 
materials  and  superlattices  provided  they  share  the  same  high  symmetry  axis 
z'.  The  validity  of  these  arguments  has  been  numerically  tested  in  I  for  a 
simple  case. 


The  second  important  ingredient  of  the  method  is  the  simplicity  of  the 
pertinent  wave  equations  for  ^^(g,z').^  This  is  a  set  of  coupled,  one¬ 
dimensional  second-order  ordinary  differential  equations  (in  z').  The 
coupling  potentials  are  derivable  from  the  underlying  three-dimensional 
pseudopotential  by  means  of  geometrical  coefficients,  the  "A-coef f icients" 

(see  I).  These  equations  are  reminiscent  of  a  set  of  coupled  Kronig-Penney 
models.^  They  can  be  solved  by  any  convenient  algorithm  such  as  straight 
integration  (when  the  energy  is  not  too  low),  or  via  conversion  into  a  secular 
matrix . 


With  spin  degrees  of  freedom  the  features  just  discussed  hold,  though 
some  of  the  details  change.  An  obvious  modification  is  that  with  spin,  the 
wave  function  is  comprised  of  two  components  pertaining  to  the  two  possible 
spin  projections.  Consequently,  two  multipole  expansions  such  as 
Equation  (1-1)  need  be  considered  simultaneously,  and  the  total  number  of 
relevant  multipoles  roughly  doubles.  The  corresponding  wave  function  for  the 
multipoles  include  additional  terms  (in  comparison  to  the  no-spin 
counterpart),  one  of  which  is  a  first  order  derivative  (Section  4).  The  SO 
effect  is  primarily  mixing  the  two  spin  components  of  the  wave  function.  In 
the  present  phenomenological  approach,  this  mixing  is  controlled  by  adjustable 
parameters.  The  relative  ease  of  including  SO  is  an  attractive  feature  of  the 
present  method. 
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To  demonstrate  how  the  method  works,  and  as  a  first  application,  we 
consider  the  PbSe  band  structure  along  the  [111]  direction.  The  analysis 
demonstrates  several  points: 

1.  Group- theoretic  considerations  directly  indicate  the  SO-induced 
admixtures.  In  particular,  we  recover  the  effect  of  band  gap  reduction  by  the 
SO  interaction. 

2.  We  correctly  obtain  the  low  bands  throughout  most  of  the  T-L  line. 
Numerical  artifacts,  however,  limit  at  present  the  method  to  non-degenerate 
bands . 


3.  The  spin-orbit  interaction  depends  on  two  phenomenological 
parameters,  roughly  corresponding  to  the  (i-0,  j-1/2)  and  (i»l,  j-1/2 ,  3/2)  SO 
strengths . 

The  paper  is  organized  as  follows.  In  Section  2  we  introduce  the 
multipole  expansion  in  the  presence  of  spin.  In  Section  3  we  discuss  the 
group- theoretic  considerations  pertaining  to  the  cylindrical  multipoles. 
Section  4  is  devoted  to  the  derivation  of  the  central  wave  equations  for  the 
multipoles.  The  last  section  describes  the  results  of  numerical  calculation 
for  PbSe  bands  in  the  P-L  direction  and  a  brief  discussion.  Most  of  the 
technical  details  are  given  in  the  appendices. 
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SECTION  2 

CYLINDRICAL  COORDINATES  MULTIPOLE  EXPANSIONS  IN  THE  PRESENCE  OF  SPIN 


The  formalism  described  in  I  is  extended  here  to  include  spin.  In  the 
unprimed  coordinates  (x,y,z),  the  spin  wave  function  is  a  spinor  (x(o))  where 
(x(l))  and  (x(-l))  designate  the  "spin-up"  and  "spin-down"  states, 
respectively.  Hereafter,  we  adopt  the  notation  that  parentheses  around  a 
function  or  operator  indicate  a  spinor,  such  as  for  (x(o)).  Using  this 
notation,  a  spinor  Bloch  function  in  the  unprimed  coordinates  (^'_T*(r))  has  the 
form^ 


(«.nij(r)) 


i  lc«T*  -► 

e  (U^.(r)) 


where  the  lattice-periodic  spinor  (Uj^]j(t))  has  the  decomposition^ 


o-±l 


(2-1) 


(2-2) 


To  derive  the  multipole  expansion  of  Equation  (2-1)  we  first  transform 
to  the  primed  coordinates  where  z'  points  at  the  direction  of  a  high -symmetry 
axis.  Such  a  rotation  affects  also  the  spinors  since 


R  ['I'(r)(x)]  =  [R  -Kr)]  (R(x)] 


(2-3) 


where  R  is  the  rotation  operator,  and  ^'(r)  and,^(x)  are  an  arbitrary  function 
and  spinor,  respectively.  The  spinor  factor  [R(x)]  transforms  according  to 
the  D^'^(a,p,'Y)  two-dimensional  matrix^’®  where  of.p.y  are  the  Euler  angles. 
The  radial  factor  [R4'(r)]  =  ^f(R'4')  has  been  discussed  in  I .  As  in  I ,  we 
limit  ourselves  to  crystal  momentum  I?'  such  that  lc'*r'  -  k'z'.  Therefore,  in 
the  primed  coordinates 


(^^^,(r')) 


ik '  z  V 

e  L 

a 


U^2>(?')(x'(o)) 


(2-4) 


where 

refer 


U^5^(r')  are  linear  combinations  of  spinors 

to  the  z ' -direction .  For  the  [111]  direction  in  particular 


(x'(o)) 
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0  <  k'  -  K  ^ 


0  <  K  <  1 


(2-5) 


The  form  Equation  (2-4)  involves  two  lattice-periodic  functions. 
Therefore  in  analogy  to  Equation  (1-1)  (see  I),  they  have  the  following 
multipole  expansion 


-I 

a=±l 


I  I  J/gpP)  (x'(o))  (2-6) 


where  the  two  spin  multipoles  are  given  by 


(2-7) 


\^\  =  gp 

In  Equation  (2-7)  denoted  a  reciprocal  lattice  vector  in  the  primed 
coordinates,  Cj  is  its  transversal  (normal  to  z')  component,  4(t')  are  the 
Fourier  coefficients  of  ^'(r'),  and  the  shift  phase  6(t.j.)  is  defined  in  I. 

We  turn  now  to  the  SO  interaction.  The  goal  is  to  cast  it  in  a  form 
compatible  with  expansions  of  the  structure  Equations  (2-6)  and  (2-7). 


In  the  unprimed  coordinates  and  cgs  units,  the  spin  orbit  interaction  is 
given  by^ ’ ^ 


/  2  2 
4mo  c 


(r)  X  V) 


(2-8) 


In  Equation  (2-8)  the  Pauli -matrices  "vector"  is  o  =  (°x’°y’°z^  V(r)  is 
the  proper  ionic  lattice -periodic  potential.  The  details  of  V(r)  are  unknown. 
However,  in  the  ionic  core  portion  of  space,  where  its  derivative  is  the 
largest,  it  is  customarily  approximated  by  a  Coulombic  potential  V^(r)  with  a 


lei 

form  factor.  For  an  ion  at  r  -  0,  it  follows  that  V  V(r)  =  V  V^(r)  -  -■ — r 

r 

where  Z  is  the  effective  ionic  charge.  Consequently,  Equation  (2-8)  takes  the 
standard  form 


r|-0  4mQ  c  r 


t _  Mfz  7) 

2  2  3 


(2-9) 
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where  the  angular  momentum  is  i=-irxv. 

For  a  lattice-periodic  array  of  ions,  a  natural  generalization  of 
Equation  (2-9)  is 


=  I  f.(l?  -  ?-!)  (o  .  i  (?  -  ?.))  (2-10) 

i 

where  r^  runs  over  all  atomic  positions  in  the  crystal  and  fj^(|t  -  ^^  |  )  is  a 
phenomenological  form  factor.  This  form  is  readily  transformed  to  the 

A  A  A  A 

primed  coordinates  since  (o  .  i)  -=  (o'  .  i')  and  distances  |r  -  r  |  are 
conserved  by  rotations.  In  addition,  since  the  nonderivative  factors  in 
Equation  (2-10)  are  lattice-periodic,  the  Fourier  expansion  of  the  SO  form 
factor  in  the  primed  coordinates  is 


f(r')  -  I  f.(|?' 
i 


-  rp  -  I  f(G') 
G' 


iG' .r ' 
e 


and  the  SO  interaction  Equation  (2-10)  takes  the  convenient  form 


(2-11) 


(Vso(?'))  -  -i  (o'  .  [  X 

G' 


f  (G*)  e 


iG'  .r 


XV']) 


(2-12) 


The  form  Equation  (2-12)  is  mathematically  convenient.  Its  operation  on 


iG  *  •  IT  ^ 

products  of  the  type  e  (x'(v)).  Equation  (2-7),  is  easily 

evaluated.  It  also  naturally  introduces  the  multipoles  of  the  SO  vector  form 
factor  Equation  (2-11),  i.e.. 


f(r')  =  I  I  f^(g  ,z')  J^(gQp) 
gQ 

where  in  analogy  to  Equation  (2-7)  (see  I,  Section  2) 


t  ,  .i. 

:/gQ.z  )  =  1 


(2-13) 


(2-14) 


The  ?^(gQ,z'),  of  dimensionality  [ E*f ] ,  are  the  SO  counterparts  of  the 
"central^  pseudopotential  multipoles  v^(gp,z'). 
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SECTION  3 

SYMMETRY  PROPERTIES 


The  assumption  that  z '  is  a  high  symmetry  axis  imposes  severe 
restrictions  on  the  allowed  i  values  in  a  multipole  expansion  such  as 
Equation  (2-6)  in  part  I.  Lattice -periodicity  in  the  transversal  direction 
(normal  to  z ' )  yields  additional  constraints  on  the  allowed  £  values.  In  this 
Section  we  consider  these  constraints  in  the  presence  of  spin  degrees  of 
freedom.  For  the  sake  of  self  containment  we  first  briefly  review  results 
from  part  I.  Technical  details  are  given  in  Appendices  A  and  B. 

Consider  the  multipoles  of  a  lattice-periodic  function  W(r')  with 
Fourier  components  W(2').,  Th^  origin  of  sjmunetry  properties  is  that  for  a 
r^cipro9al  yector  2'  -  (G^,  Gy,  G^)  in  the  primed  coordinates,  all 
^T  “  ^^x’  ^y^  such  that  |G^|  =  g,  lie  on  spokes  of  one  or  two  conjugated 
hexagons  (Figure  3-1,  I).  Furthermore,  there  are  simple  phase  relations 
between  W(2')  and  W((2')j^)  where  (2')f^  -  ((2j)j^,  G^ ' ) ,  (^x^r  -  (G^'.  -  G  '). 
These  properties  and  the  structure  of  a  multipole  Equation  (2-7)  lead  to'^a 
situation  in  which,  for  a  particular  £  value,  either  all  terms  in  the 
summation  Equation  (2-7)  cancel  out  or  add  up.  As  a  specific  example,  the 
pseudopotential  v(r),  by  its  very  construction,  is  invariant  under  all 
symmetry  operations  of  the  fee  lattice.  In  particular,  it  transforms  as  the 
scalar  under  the  0^  group  and  under  all  subgroups  of  0^,  such  as  the  C^.^  and 
Djjj  groups  associated  with  the  [111]  direction. This,  in  turn,  restricts 

the  i-values  in  Equation  (1-1)  to  i  -  3m,  m  —  0,  +1 . (I,  Table  3-3).  In 

addition,  by  virtue  of  the  structure  of  the  multipole  v^(gp,z'),  it  is  £- 
periodic  ( see  I ) : 


^£±6  v^(gp,z') 

for  sin  [6$Q(g)]  -  0  (3-1) 

where  $Q(g)  is  given  in  Table  1-1.  These  two  constraints  imply  that  i  -  0,  3 
are  the  only  independent  multipoles  for  v(r);  all  others  are  phase  related. 

In  the  presence  of  spin,  the  underlying  symmetry  groups  are  the  double 
groups  and  (see  Appendix  A).  The  corresponding  £  sequences  are 

derived  in  Appendix  B  and  summarized  in  Tables  3-1,  3-2.  Note  that  the 
allowed  £  values  for  the  "up"  and  "down"  spin  components  are  different .  albeit 
they  still  belong  to  sequences  with  an 
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TABLE  3-1.  MULTIPOLE  CONTENT  OF  THE  IRREDUCIBLE  SPIN  REPRESENTATIONS  OF  03^ 

i -parity 

IRREDUCIBLE 

REPRESENTATION  i -values ^ 


A4 

1 1>“ 

■3m  -11 
.3m  +lj 

(.l)3m 

0 

-i 

i' 

0. 

1 

^5 

|1>=  1 

'3m  -1] 

.3m  +lj 

.(.l)3m 

0 

-i 

i' 

0. 

1 

^6 

ii>- 

(-l)3m 

0 

0 

0 

0 

0 

1 

1‘ 

0 

2 

0 

1 

0 

0 

1 

0 

0 

0 

(a)  m  -  0,  +1,  +2,  +3,  .... 


(b) 

'  U_^(g,z') 

A 

=  ‘Tl  « 

■  u^(g.z')  ' 

V_^(g,z') 

i 

V/g,z') 

■  U_^(g.z')  ■ 

■  U^(g.z') 

or 

V_^(g,z') 

A 

v^Cg.z') 

W_^(g,z') 

W/g.z') 

X.^(g.z') 

X^(g.z’) 

NOTE:  The  i- sequences  pertaining  to  the  representations  of  the  double  group 
C^J) .  The  notation  is  explained  in  Appendix  B. 
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TABLE  3-2. 


MULTIPOLE  CONTENT  OF  THE  IRREDUCIBLE  SPIN  REPRESENTATIONS 


OF  D 


(d) 

3d 


IRREDUCIBLE 

REPRESENTATION  i -values 


i-parity  z-parity 


d 


|i>= 

'3m 

3m 

(.l)3m 

’0  -i 

.i  0. 

(-1)^ 

ri  o-i 
0  ij 

|i>= 

3m 

3m 

0  -i' 
i  0 

1  01 
0  ij 

|i>= 

3m 

3m 

-r 

+1, 

(-l)3m 

’0  -1 
i  0 

-(-1)^ 

1  01 
0  ij 

|i>= 

3m 

3m 

-r 

+1. 

-(-l)3m 

0  -i 
i  0 

-(-1)^ 

1  01 
0  ij 

|1>= 


3m 

,3m  +1, 


1,3m  +lj 


■3m  -r 

'0 

0 

0 

1' 

'1 

0 

0 

O' 

(.l)3m 

0 

0 

1 

0 

0 

-1 

0 

0 

3m 

0 

1 

0 

0 

0 

0 

-1 

0 

,1 

0 

0 

0. 

,0 

0 

0 

1 

'0 

0 

0 

r 

’1 

0 

0 

O' 

3m  -1 

(.l)3m 

0 

0 

1 

0 

-(-l)3m 

0 

-1 

0 

0 

0 

1 

0 

0 

0 

0 

-1 

0 

3m  J 

1 

0 

0 

0. 

,0 

0 

0 

1, 

2 


2 


(a)  m  =  0 ,  +3 ,  +6 ,  .... 

A 

(b)  See  Table  3-3  for  the  definition  of  n^. 


(c) 

■  U/g,-z')  ■ 

A 

■  U^(g,z') 

■  U^(g.'Z')  ■ 

■  U^(g,z') 

V^(g. -z') 

S  It 

z 

V^(g.z') 

or 

V^(g.-z') 

W/g.-z') 

X/g.-z') 

> 

N 

V^(g.z') 

W^(g.z') 

NOTE:  The  i-sequences  pertaining  to  the  representations  of  the  double  group  The 

notation  is  explained  in  Appendix  B. 
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increment  of  ±3.  Equation  Equation  (3-1)  is  valid  for  each  spinor  component 
individually.  When  sin[ 6$Q(g) ]?'0  the  i -periodicity  relationship  becomes  more 
complex,  however,  this  situation  does  not  arise  in  our  example  (Section  5 
and  I ) . 


The  comparison  of  the  multipole  content  of  the  wave  function  (i.e.,  the 
allowed  f  values)  with  spin  (Tables  3-1,  3-2)  with  the  corresponding  tables  in 
the  no-spin  case  (Tables  3-2,  3-3  in  I)  is  instructive.  By  matching  the 
allowed  i-value  sequences,  it  is  possible  to  infer  the  representations  admixed 
by  the  SO  interaction.  Consider  first  a  A-point,  i.e..  Table  4-1  and 
Table  3-2  in  I.  From  the  comparison  we  conclude  that  the  spinless  A^,  A2 
representations  are  admixed  with  the  A^  representation  to  yield  the  A^  spin 
representation.  On  the  other  hand,  the  spin  representations  A^,  A^  result 
from  the  splitting  the  degeneracy  of  the  spinless  A^  representation  by  the  SO 
interaction.  At  the  L-point,  comparing  Table  3-2  with  Table  3-3  in  I  implies 
that  the  SO  interaction  mixes  ,  Ln,  with  (and  '  ,  L2  '  ,  with  ' )  to  give 
the  spin  representation  Lg'*'  (or  L^'). 

Consider  now  the  SO  vector  form  factor  f(r').  Equation  (2-11).  By 
comparing  it  to  Equation  (2-8),  it  follows  that  f(r')  -  ^'V(r')  where  V(r') 
transforms  as  a  scaler  under  the  point  group  transformations.  Since  V(r')  has 
the  multipole  expansion  Equation  (1-1)  with  i  -  3m,  m  -  0,  ±1,  ±2,...,  and 


d 

-  COS<|) 

sin4> 

a 

ax' 

dp 

P 

a<)> 

a 

-  sin4) 

a  ^ 

COS()> 

P 

a_ 

ay' 

dp  * 

8^ 

a 

a 

az'  az' 

(3-3) 


it  follows  from  the  Bessel  functions  properties^  that  the  only  SO  non¬ 
vanishing  components  are 

^L+1  ^  ^3m+i  ^3m  ^ 

^3m-l  ^  ^  ^3^1  ^  ^3m  ^ 


(gq.z  ^  “  dz'  ^3m  ^ 


■3m 


m  -  0 ,  ±1 , 


±2, 


(3-4) 


where  =  (f^^  \gQ,z'), 

are  the  multipoles  of  V(r) . 


f^y'^gp.z'),  f^^'\gQ,z'))  and  V3jgQ,z') 
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Eliminating  '^3in(gQiZ')  from  Equation  (3-4)  yields  general 
relations  for  the  SO  vector  form  factor: 

f^^  ^  (g  z')  =  -  f^^  ^  (g  z')  -  i  f^^  ^  (g  z') 

3m+l  ^  3m-l  ^  3m+l  ^ 

f^y  ^  (g  z')  -  f^^  ^  (g  z')  -  i  f^^  ^  (g  z') 
^3m+l  '‘®Q’  ^  ^3m-l  ''  3m-l  '^^Q’  '' 


When  the  i-se 
follows  that 


for  m  -  0 ,  ±1 ,  ±2 ,  ... 

que^ce  of  Equation  (3-5)  is  compared  to  Table  3-3 
f(x  f(y  ^(z')  transform  as  Lo '  while  f^^  ^( 


as  the  1,2'  representations  of  D^^j. 


symmetry 


(3-5) 

in  I ,  it 
z')  transforms 
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SECTION  4 

THE  MULTI POLES  WAVE  EQUATIONS 


In  this  section  the  wave  equation  for  the  multipoles^  is  generalized  to 
include  the  spin  degrees  of  freedom  and  the  SO  interaction.  We  follow  the 
method  in  I ,  i.e.,  first  the  interaction  terms  are  transcribed  to  a  product 
form  and  then  the  decomposition  theorem  is  applied. 

In  the  primed  coordinates,  the  Schrodinger  equation  with  spin-orbit 
interaction  is 


[-V'2  +  E  (k) 


2m_  2m-  ^ 

V(?')  -  (VgQ(?'))](4r^-(?'))  =  0 

h  h 


(4-1) 


A 

where  (VgQ(r'))  is  the  SO  interaction  and  mQ  is  the  effective  mass,  chosen  to 
be  the  free  electron  mass.  Inserting  the  wave  function  multipole  expansion 
Equation  (2-6)  into  Equation  (4-1),  projecting  on  e^**  and  factoring  out  the 
Bessel  functions,  we  obtain: 


I 

o-±l 


2ni 

+  E^(k)  + 


(o) 

i 


(g 


2m 

o 


(gQ-Z  ) 


2m 


(x'(o))  “  0 


(4-2) 


In  Equation  (4-2)  t^je  and  (Y^^°^)  terms  originate  from  the  products 

"v(r '  )'I'(r ' ) "  and  "(V5Q(r')  'ff(r'))"  in  the  Schrodinger  Equation,  respectively. 
The  former  has  been  discussed  extensively  in  I.  The  latter  is  new. 
Specifically,  the  goal  is  to  express  (Y^'°^ (gQ, z ' ) )  in  terms  of  multipoles 
pertaining  to  the  SO  form  factor  and  the  wave  function.  Details  are  given  in 
Appendix  C;  the  rest  of  the  details  are  readily  carried  over  from  I.  These 
steps  yields  the  multipoles  wave  equation. 
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fi  dz 


+ 1 


3 

I 


gp-Sp 


[  ^"*0 

Vm(gp,z')  0- 

« - * 

f 

»  1 
N> 

0 

2m, 


0 


dz'  i-m'‘®F’  ’ 

z') 

az'  i-m'‘®F’  ' 


^"*0  ^F 

^7 


^"*0  £f 

ZT  ~ 


■ 

0 

f  (R  z')  1 

•2f‘^'>(gp..') 

.  ’*'i-m-l^®F’^  ^ 

*i-m+l^®F’^  ^ 

0 

(4-3a) 


The  symbols  in  Equation  (4-3)  denote  the  following:  The  band  energy  is 
denoted  by  Ej^(k'),  the  (gp.z ' )  are  the  two  spin  components  multipoles  of 

■‘  ')),  ■v„(gp.z')  are  the  multipoles  of  Y(r')  - 
(4-1)),  the  SO  vector  multipoles  f^“'(gp,z') 


the  wave  function  (Equation  (2-6 
the  central  potential  (Equation  ,  . 

and  f^^  (gp.z')  are  defined  in  Equations  (2-14),  (3-4)  and 


(4-3b) 


and  the  geometrical  ^§0’  gp’  gp) 'Coefficients  have  been  introduced  and 

discussed  in  I.  The  selection  rules  for  i,  m,  gp  and  gp  are  incorporated  into 
the  A- coefficients .  The  boundary  conditions  for  the  wave  function  multipoles 
are  derived  from  the  Bloch  theorem  in  the  [111]  direction 

t(o),  ,  ik'a*,  /  //  1  \ 

^(gp,  z'  +  a  )  -  e  (gp,  z')  (4-3c) 

where  k',a  are  defined  in  Equation  (2-5). 

The  single  variable  wave  equations.  Equation  (4-3),  is  the  central 
result  of  this  work.  The  band  symmetry  is  specified  by  the  selected  i- 
sequence  according  to  Tables  3-1  and  3-2.  Note  that  the  i- sequence  for 
4'^°^(gF,z')  depends  on  o,  the  spin  projection.  The  selected  g  values,  in 
ascending  order  (Table  1-1),  are  controlled  by  considerations  of  convergence 
and  available  computing  capacity. 
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An  interesting  feature  of  Equation  (4-3)  is  the  appearance  of  first 
derivative  (3/3z')  terms.  They  originate  from  the  z'  derivative  terms  in  J.s. 
The  SO  interaction  in  Equation  (4-5)  has  both  diagonal  and  off-diagonal  terms 
in  o.  The  former  snift  the  no-spin  bands  as  determined  by  the  central 
potential  ■V[„(gp.z')  while  the  latter  admix  or  split  the  no-spin  bands.  While 
in  Equation  (4-5)  both  terms  are  controlled  by  the  same  spin  interaction,  we 
will  argue  in  the  next  section  that  in  the  context  of  a  phenomenological 
approach  these  two  types  of  terms  should  be  parametrized  separately. 
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SECTION  5 

BAND  STRUCTURE  OF  PbSe  IN  THE  F-L  ([111])  DIRECTION 


To  demonstrate  the  viability  of  the  present  method,  we  calculate  in  this 
section  the  low  bands  of  PbSe  along  the  F-L  direction.  This  example  has  been 
analyzed  in  I  in  the  absence  of  SO  interaction.  The  analysis  serves  a  dual 
purpose:  to  examine  the  implementation  of  the  wave  equation  (4-3)  and  to 
demonstrate  the  type  of  physical  insight  the  present  approach  provides.  Since 
the  convergence  of  the  multipole  expansion  has  been  demonstrated  in  I,  we 
consider  here  only  the  SO  effects.  To  keep  the  discussion  focused,  the 
section  is  subdivided  according  to  the  various  aspects  of  the  analysis. 

THE  SPIN-ORBIT  INTERACTION 

We  discuss  first  Equation  (4-5)  in  the  context  of  the  Lg(Ag)  and 
^4  5(4^  5)  symmetries.  The  i -values  for  which  the  SO  form  factor 
Equation  (2-14)  does  not  vanish,  given  the  symmetry  relations  Equation  (3-5), 
are  shown  in  Table  5-1.  The  wave-function  i-sequences  for  the  above 
symmetries  are  given  in  Tables  3-1,  3-2.  Combining  these  two  i-constraints 
yields  the  wave  equation  block  structure  in  o,  the  spin  projection,  displayed 
in  Table  5-2:  For  the  Lg(Ag)  bands,  the  SO  terms  occur  in  one  diagonal  block 
(i  =  3m+l)  and  the  off  diagonal  blocks,  Table  (5-2a).  On  the  other  hand,  for 
the  L2|  5(4^  5)  bands  only  the  diagonal  blocks  (i  -  3m±l)  have  SO  terms. 

Table  (5-2b).  This  structure  has  a  physical  interpretation  recalling  that  the 
i  -  3m  blocks  are  associated  with  s-waves  (and  higher  even  waves)  whereas  the 
St  =  3m±l  blocks  with  p-waves  (and  higher  odd  waves).  Thus  the  diagonal  SO 
terms  in  the  St  =  3m±l  describe  the  p-wave  SO  interaction  while  the  off 

diagonal  SO  terms  describe  the  s-p  SO  mixing.  Since  the  ^1/2  Pl/2 

interactions  are,  in  general,  quite  different, it  is  expected  that  these  two 
types  of  terms  have  different  strengths. 

The  parametrization  of  the  SO  interaction  depends  on  the  specifics  of 
the  material.  Since  PbSe  has  two  atoms  in  the  unit  cell,  the  SO  form  factor 
f.j(|r-r  J),  Equation  (2-11)  and  Appendix  D,  requires  four  parameters:  the  two 
strengths  and  two  ranges  of  the  atomic  SO  interactions.  The  SO  ranges  are 
taken  to  be  that  of  Pb  and  Se  atomic  cores;  The  Pb  core  (Xe  4f^^  5d^®)  is 
approximately  that  of  Ar  and  the  Se  core  (Au  3d  )  is  approximately  that  of 
Cu.  The  core  radii  values  are  given  in  Table  5-3.  The  radial  dependence  is 

taken  to  be  a  gaussian.  A  box  form  factor  gives  similar  results  (Appendix  D) . 
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TABLE  5-1.  SYMMETRY  PROPERTIES  OF  THE  SPIN-ORBIT  FORM  FACTOR  f,n(gp.z') 


(z')  (b) 


(a)  A(-z')  -A(z'),  B(-z')  -  -B(z') 

(b)  C(-z')  -  -C(z'),  D(-z')  -  D(z') 

NOTE:  The  non-vanishing  components  and  symmetries  of  the  SO  form  factor 
(Equation  (3-4),  (4-3b)). 
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TABLE  S.2.  BLOCK  STRUCTURE  OF  WAVE  EQUATIONS  (4-3) 


0--  1 

0-+1 

(3m+l) 

(3m) 

o=-l 

(3m+l) 

v„(gp,z')  +  fl*'(gp,z') 

o=+l 

(3m) 

v^(gp,z') 

(a);  Lg(Ag)  representations 


0--1  0-+1 

(3m+l)  (3m-l) 


0--  1 
(3m+l) 


Vn,(gp.z')  +  f^-\gp,z') 


0 


0-+1 

(3m-l) 


0 


(b):  L^  5(^4  5)  representations 


NOTE:  Schematic  block  structure  of  the  wave  functions  Equation  (4-3)  for  the  L 

and  L2|  5(^4  5)  bands.  The  o  and  i-sequence  indices  are  explicitly  indicated. 
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TABLE  5-3.  SPIN-ORBIT  PARAMETERS 


SO 


s(  o  •  i  ) 


Ase  [eV]  Apj^  [eV]  r(Pb)(A)  r(Se)(A) 


REFERENCE 


J.  C.  Philips^^) 

.16 

.66 

H.  J.  Zeiger 
and  G.  W.  Pratt 

.14 

.45 

This  Work 

.14 

.56 

1.45 

1.16 

(a)  We  use  A  -  .33*A^  where  is  a  spin-orbit  splitting  of  Atom  A, 

Equation  (7-28)  in  Reference  3,  p.l79. 

(b)  We  use  Equation  (3-61)  in  Reference  7  for  the  Atomic  Spin-orbit  Strengths. 


NOTE:  The  spin  orbit  parameters,  see  Appendix  D.  The  s  and  r  parameters  denote  the 
atomic  SO  strength  and  range  parameters,  respectively. 


5-4 


NAVSWC  TR  91-326 


The  SO  strengths  are  chosen  as  follows.  The  ratio  of  the  Pb  and  Se  SO 
strengths  is  taken  from  atomic  data  (see  Table  5-3).  This  ratio  is  roughly 
3.-4.  and  is  kept  fixed.  The  overall  SO  strengths,  however,  are  adjustable 
parameters.  According  to  the  above  remark,  the  diagonal  SO  blocks  (in  Table 
5-2)  pertain  to  the  p-wave  whereas  the  off  diagonal  SO  blocks  pertain  to  s-p 
SO  mixing.  Consequently  we  introduce  two  dimensionless  SO  scale  parameters: 
SjjQ,  a  nondiagonal  SO  parameter  and  Sjj  an  diagonal  SO  parameter  such  that 


A(Pb)  = 


Pb’ 


A(Se) 


Se 


(5-1) 


where  Ap^,  are  the  "atomic"  spin  orbit  strengths  of  Pb  and  Se ,  Table  5-3, 

and  A(Pb),  A(Se)  are  the  overall  atomic  SO  strengths  (see  Appendix  D) .  The 
signs  and  magnitudes  of  S^  are  determined  by  fitting  two  measured  band 

gaps  (see  the  next  subsection).  Consider  the  Lg  bands  (Table  (5-2a).  The  "v" 
terms  in  the  i=3m  block  correspond  to  the  s-^^2  potential.  By  the  same  token, 
the  "v  +  f“"  terras  in  the  ^=3ra+l  block  correspond  to  the  potential. 

Since  the  latter  is  less  positive  than  the  former, we  can  expect  Sjj<0  to 
make  "f-"  negative.  Strictly  speaking,  the  wave  equation,  Equation  (4-3), 
implies  that  Sj^j^  =  Sj^.  However  if  it  is  reasonable  to  expect  them  to 

share  the  same  sign.  Thus,  heuristically ,  we  expect  both  Sj^j^  and  Sj)  are 
negative . 


RESULTS  FOR  PbSe  ALONG  THE  [111]  DIRECTION 

The  four  relevant  bands  near  the  Fermi  level  are  apparent  from 
Figure  5-la.  These  particular  no-spin  bands  are  calculated  using  an  eight- 
multipole  basis,  which  is  a  computational  compromise  between  a  fully 
convergent  basis  (12  multipoles)  and  our  computing  capacity.  A  consequence  of 
the  eight-multipoles  basis  choice  is  that  the  no-spin  bands  are  not  adequately 
described  for  K<.5  (Equation  (2-5)),  i.e.,  for  the  lower  half  of  the  F-L  line. 

t  r 

At  the  L- point,  the  SO  interaction  admixes  the  (L-j,L]^)  and  the  (L2,L2) 
bands,  respectively.  This  is  evident  from  comparing  the  ^-sequences  for  the 
Lg  bands  (Table  3-2)  with  those  of  L^  and  Lj^  bands, (see  I).  Consequently, 
given  the  particular  arrangement  of  the  ,L-j  ,L2^ , L2 , L^  bands  in  Figure  5-la,  it 
is  expected  that  the  band  gap  AE^  =  E(L2)  -  E(L^)  will  be  reduced  by  the  SO 
interaction.  This  trend  is  demonstrated  in  Table  5-4  and  Figure  5-2  for  the 
near-Ep  Lg  bands.  Consider  first  the  limit  of  Sjjp  -  0.  From  Table  5-2  it 
follows  that  in  this  case  only  the  p-wave  (£-3m+l)  blocks  are  modified.  This 
is  evident  from  the  first  two  columns  of  Ta^le  5-4  where  the  L^.L^  bands  (the 
two  outermost  Lg  bands)  move  while  the  Lp,L2  bands  (the  two  innermost  Lg 
bands)  do  not.  Positive  Sq  tends  to  cluster  the  bands  while,  conversely,  Sq<0 
fans  out  the  bands.  In  the  other  limit,  i.e.,  Sjj  -  0,  Table  5 -2a  implies  that 
only  off-diagonal  s-p  mixing  blocks  exist.  In  this  case  the  band  gap  is 
reduced  as  expected  and  all  bands  move  down.  This  statement  is  true  for  both 
positive  and  negative  Uhen  both  and  Sjj  are  nonzero,  a  combination  of 

the  above  two  effects  is  observed.  In  particular,  when  and  have  the 
same  sign,  then  the  bands  either  cluster  together  for  or  fan  out 
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TABLE  5-4.  BAND  ENERGY  DEPENDENCE  ON  SPIN-ORBIT  PARAMETERS  AT  THE  L-POINT  Lg 


Band 


^ND 

Sd 

Symmetry 

0 

0 

0 

-20. 

0 

20 

±30. 

0 

30. 

30. 

-30. 

-30. 

^6 

10.718 

11.436  t 

9.965  1 

10.573  t 

9.977  1 

11.670  t 

9.424 

9.424 

9.424 

9.312  1 

9.439  t 

9.297  i 

^6 

8.853 

8.853 

8.853 

8.830  1 

9.093  1 

8.792  1 

^6 

7.180 

6.816  1 

7.503  t 

7.053  1 

7.501  t 

6.540  1 

AE 

0.571 

0.571 

0.571 

0.482 

0.366 

.505 

NOTE:  The  dependence  of  near-Ep  PbSe  band  energies,  at  the  L-point,  on  the  two  spin- 
orbit  parameters,  Equation  (5-1).  The  trend  of  the  band  movement  (in  energy) 
is  indicated  by  an  arrow  and  the  band  gap,  AE.  All  energies  are  given  in  eV. 
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for  Sjj£j,Sq<0  and  the  band  gap  is  reduced.  Figures  5-2a  and  5-2b  demonstrate 
the  trend  for  the  physically  relevant  case  when  Note  that  a 

sufficiently  strong  SO  interaction  reduces  the  band  gap  to  the  point  of  "band 
inversion,"  and  beyond.  This  effect  of  the  SO  interaction  is  well  known  in 
narrow  gap  semiconductors,  such  as  in  the  IV-VI . 

Figures  5 -lb  and  5-lc  show  the  calculated  lowest  ten  bands  and  the 
corresponding  results  from  Kohn's  work.^^  Numerical  details  of  the 
calculation  are  given  in  Appendix  E.  We  choose  the  SO  parameters  as  Sj^p  = 

-60.  ,  Sjj  ■=  -15  to  approximately  reproduce  two  experimental  band  gaps  at  the  L- 
point:  the  smallest  gap  (E(L^)  -  E(Lg))  and  a  secondary  gap  between  the  two 

adjacent  Lg  bands.  As  Figures  5 -lb  and  5-lc  show,  this  choice  reproduces 
quite  well  all  the  bands  and  relative  spacings  at  the  L-point  and  down  to 
about  the  raid-  point  of  the  F-L  line.  For  smaller  K  values,  some  bands  are 
missing,  and  not  all  bands  are  in  the  correct  order  (Figure  5-lc  is  the 
correct  band  structure).  We  understand  these  deficiencies  as  numerical 
artifacts  that  originate  from  the  eight-multipoles  no-spin  base  choice.  As 
Figure  5-la  shows,  for  K  =  .  5  we  have  a  spurious  crossing  of  the  and 

bands.  This  crossing  creates  a  wrong  sequence  for  two  bands  at  the  F-point. 
Obviously,  the  SO  mixing  cannot  rectify  this  deficiency;  hence,  the  wrong 
sequence  at  the  F-point  (Figure  5-lb).  In  addition,  the  (Lj^,  L^)  crossing 
generates  a  degeneracy  of  equal -parity  bands.  For  numerical  reasons,  the  code 
does  not  deal  correctly  with  a  situation  of  several,  closely  spaced  and  mixed 
bands.  This  is  the  reason  for  the  missing  of  a  few  bands  for  limited  range  of 
K  values.  The  same  situation  occurs  at  E~14eV,  where  again  two  closely  spaced 
bands  occur. 

From  the  comparison  of  Figures  (5-la)  through  (5-lc),  we  conclude  that, 
whenever  the  no-spin  basis  is  adequate  (i.e.,  K  «  .5  in  our  case),  the  bands 
calculated  by  the  present  method  compare  well  with  the  correct  band  structure. 
Conversely,  whenever  the  no-spin  basis  is  deficient  (i.e.,  K  . 5  in  our 
case),  the  calculated  bands  with  the  SO  interaction  are  deficient.  This 
outcome  obviously  suggests  the  choice  of  the  fully  convergent  no -spin  base  as 
a  starting  point,  i.e.,  the  12  multipoles^  (valid  throughout  the  F-L  line) 
rather  than  the  eight-multipoles  basis  used  here.  Such  a  calculation, 
however,  is  beyond  our  present  numerical  capacity. 
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SECTION  6 

SUMMARY  AND  DISCUSSION 


We  have  generalized  the  cylindrical  multipole  approach  to  band  structure 
to  include  spin  degrees  of  freedom  and  the  SO  interaction.  The  ensuing  one¬ 
dimensional  set  of  wave  equations  is  solved  for  the  test  case  of  PbSe  in  the 
[111]  direction.  The  results  ar|^compared  with  the  band  structure  calculated 
by  an  entirely  different  method. 

The  combination  of  the  formalism  and  test  calculations  demonstrates 
advantages  and  limitations  of  the  present  approach.  We  have  shown  that  the  SO 
interaction  is  accommodated  naturally  in  the  present  formalism.  Another 
advantage  of  the  present  representation  is  that  it  obviates  the  band 
symmetries  which  are  admixed  by  the  SO  interaction.  Thus,  for  instance,  we 
can  readily  understand  the  band-gap  reduction  at  the  L-point  by  the  SO 
interaction  and  the  tendency  to  band  inversion  for  sufficient  strong  SO 
interaction.  The  calculated  PbSe  band  structure  in  the  [111]  direction  yields 
the  correct  sequence  and  spacings  at  and  around  the  L-point,  where  the  Fermi- 
level  lies. 

On  the  other  hand,  the  method  in  its  present  form  does  not  reproduce, 
degenerate,  or  quasi -degenerate  bands.  At  the  F-point  this  problem  can  be 
possibly  alleviated  provided  the  F-point  symmetries  can  be  introduced.  At 
other  points  on  the  F-L  line,  a  fully  convergent  no-spin  basis  avoid  incorrect 
band  ordering. 
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APPENDIX  A 

IRREDUCIBLE  REPRESENTATIONS  OF  THE  DOUBLE  GROUPS  Cj^W) 

Double  groups  arise  when  spatial  symmetry  operators  act  on  products  of 
the  form  ^(r)(x)  where  (y)  is  a  spinor.  Since  (Equation  (2-3) 

A  A  A 

R  [’i'(r)(x)]  -  [R  'l'(r)]  lR(x)]  (A-1) 

A 

the  matrices  of  R  in  this  basis  form  a  group,  termed  the  "double  group." 

A 

The  term  reflects  the  fact  that  [R(x)]  changes  sign  when  any  Euler  angle  is 

A 

changed  by  2n.  Hence  to  every  rotation  R  correspond  two  matrices  (R)  and  (R) 

such  that  (R)  -  Thus  the  order  of  a  double  group  is  twice  that  of  the 

original  group. ^ 

For  the  sake  of  completeness  we  quote  again^’^  the  cylindrical 
coordinates  realization  of  the  groups  C^.^.  and  D3JJ: 


A-1  Morgan,  D.  J.,  Solid  State  Theory.  Landsberg,  P.  T. ,  Ed.,  Wiley- 
Interscience ,  London  1965,  pp.  232-254. 

A- 2  Slater,  J.  C. ,  Quantum  Theory  of  Molecules  and  Solids.  Vol .  2,  McGraw- 
Hill,  New  York,  NY,  1965,  Appendix  9. 

A- 3  Agassi,  D.  and  Restorff,  J.  B. ,  Pseudopotential  Band  Calculations  Along 
a  Hizh- Symmetry  Axis:  Part  I- -Central  Potential  and  the  [111] 
Direction.  NAVSWC  TR  91-324,  15  Jun  1991). 
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^3v 


'3d 


Xo'I'(<|))  - 

-  'I'((|)±2ti/3) 
Yo'P(4>)  =  'K-t) 
Y^l'I'(4.)  -  'K-4.±2Tt/3) 


Xo'I'(4».z)  -  *(((..2) 

X*j^*(4),z)  -  'l'(4>±ii/3  , -z) 

Xi:24'(4.,z)  -  4'(4>±2ti/3,z) 

X34'(4»,2)  »  ^’((p+n, -z) 

Yo'I'(<|).z)  -  *(-<(>. z) 

Y^lf(<f>.z)  -  1'(-(|>±it/3,-z) 

Yi2«'(4>.z)  -  fl'(-4>^2ii/3,z) 

Y3'K(}..z)  -  'I'C-ffr+Ti.-z)  (A-2) 

In  Equation  (A-2)  is  an  arbitrary  function  and  4>.z  are  the  cylindrical 
azimuthal  angle  and  z  coordinate,  respectively.  Note  that  Yq  is  a  plane 
reflection  and  X3  is  the  inversion  operator  with  respect  to  the  origin. 

According  to  Equation  (A-1),  to  calcylate  the  double  group  irreducible 
representations  requires  the  knowledge  of  R(x)-  This  matter  is  discussed  in 
Reference  (A-2)  and  we  bring  here  only  the  results.  For  the  double 

group  corresponding  to  €3.^,  the  results  are: 


,(d) 

'3v 


>  -  n"^|a  > 


X,i|a 

Yg  |a  >  =  |p> 

Y 


±1 


ill 

a  >  -  n  o(  > 


x,i|p>  - 

Yg  |p>  -  -la  > 


^illP> 


n  -  e 


iTi/3 


(A-3) 


where  in  Equation  (A-3)  and  throughout  this  appendix,  (a)  and  (p)  denote  "spin 
up"  and  "spin  down"  spinors,  respectively. 

A 

The  representations  F  are  defined  according  to 


R  |i>-Z  |j><j  |R|i  > 
j 


(A-4) 


,(d) 


where  |  i>  is  an  appropriate  basis  state.  For  the  irreducible  spin 

representations,  we  choose  the  following  basis; 


A-2 


NAVSUC  TR  91-326 


■4- 

|1  >  -  e^'^'la  >  +  i  e‘^'**|p  > 

5- 

|l  >  -  e^*^|a  >  -  i  e'^‘^|p  > 

6- 

A 

CO. 

II 

A 

CN 

A 

II 

A 

(A-5) 

The  ensuing  irreducible  representations  are  given  in  Table  A-1. 

The  group  03^^  the  double  group  corresponding  to  contains  in 

addition  to  the  operators  Equation  (A-3)  the  space  inversion  and  a  ±n/i 
rotation  of  4).  The  former  does  not  effect  spinors.  Using  the  techniques  of 
Reference  A-2,  we  obtain: 


^.1 

|a  >  - 

|a  >  , 

x,i|p>  - 

n-^|p> 

X3 

|a  >  - 

|a  > 

X3  |p>  - 

|P> 

ja  >  ■ 

-  >  , 

1P>  -  -n"^|a> 

^0  1 

|a  >  = 

IP  > 

0 

V 

-  |a> 

(A-6) 


The  irreducible  representations  in  Table  A-2  are  obtained  from  Equation  (A-6) 
and  the  following  basis: 


Lp 


w- 


L5: 

Lt: 


1  >  -  (e^*^|o(  >  -  i  >]  U(z) 

1  >  -  (e^‘^|a  >  +  i  e‘^‘^|p  >]  U(2) 

1  >  =  [e^'<’|a  >  -  i  e’^'^lp  >]  U(z) 

1  >  -  [e^‘^|a  >  +  i  e*^‘^|p  >]  U(z) 

1  >  -  [a  >  U(z)  ,  |2  >  -  |p  >  U(z) 

1  >  -  |a  >  U(z)  ,  |2  >  -  |p  >  U(z) 


U(z)  -  -U(-z) 
U(z)  -  -U(-z) 
U(z)  -  U(-z) 
U(z)  -  U(-z) 
U(z)  -  U(-z) 
U(z)  -  -U(-z) 


(A-7) 


The  basis  functions  Equation  (A-5)  and  (A-7)  manifest 
projection  in  the  z'  direction,  denoted  by  mj  in  Tables  A-1 


the  spin- 
and  A-2. 
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TABLE  A-1. 


IRREDUCIBLE  SPIN  REPRESENTATIONS  OF  C 


(d) 

3v 


REPRESENTATION:  A^  A3 


(a) 


mj:  ±3/2  ±3/2 

OPERATOR 

Xq  1  1 

-1  -1 

Yq  -i  i 

(Plane  reflection) 


Y 


i 


±1/2 


(a)  n  - 


NOTE:  The  irreducible  spin  representations  of  the  double  group  C3v  •  The 

group  operators  are  defined  in  Appendix  A.  The  designation  of  the  spin 
projection  in  the  z-direction,  mj,  is  predicated  on  the  corresponding 
basis  functions,  see  Appendix  A. 
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TABLE  A-2.  IRREDUCIBLE  SPIN  REPRESENTATIONS  OF  D 

id 


REPRESENTATION:  L^ 

^5 

q  L5 

L^  (a) 

Lg 

mj:  ±3/2 

±3/2 

±3/2  ±3/2 

±1/2 

±1/2 

OPERATOR 


Xq  1 

111 

[J  ?]  [i 

?] 

x,i  -1 

-111 

O’)  r-n*^ 

[  0  [  0 

X±2 

-1  -1  -1 

O'! 

[  0  [  0 

01 

X3  1 

1  -1  -1 

ri  0]  r-i 

(z  parity) 

[0  ij  [0 

-ij 

Yq  i 

-i  i  -i 

ro  -i]  ro 

(Plane  reflection) 

U  oj  Li 

oJ 

Y±1  -i 

i  i  -i 

r  5,  -n'^1  r  0 

0  J 

Y±2  -i 

i  -i  i 

r  0  r  0 

[n-i  0  J  [n*i 

-Q*h 

0  J 

Y3  i 

-i  -i  i 

ro  -1]  r  0 

U  oJ  [-1 

oJ 

X.  .  -Xi  .  Yi  -  -Yi 

(a)  n  - 

NOTE:  The  irreducible 

representations  of 

the  double  group 

The  group 

operators  are  defined  in  Appendix  A.  The  mj  symbol  is  the  same  as  in 
Table  3-1. 
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APPENDIX  B 

MULTIPOLE  CONTENT  OF  THE  AND  SPIN  REPRESENTATIONS 


The  allowed  i.  values  for  the  spin  representations  of  and  the  D^^^ 

groups  are  determined  by  the  method  used  in  Reference  B-1.  To  demonstrate  the 

method,  we  analyze  here  the  Ag  representation  of  the  C^^^group.  A  similar 
analysis  applies  to  all  other  representations.  The  results  are  given  in 
Tables  3-1  and  3-2  of  the  main  text. 


two . 


From  Table  3-1  it  follows  that  the  Ag  representation  is  of  dimension 
The  general  multipole  expansion  of  the  two  basis  spinors  is 


|l>  “I  Uo  (P.z) 


,ii<l) 


|a  >  +  X  V»  (p,z)  e 
i--® 


iiij) 


Ip  > 


|2>  -  X  (p,2)  e^-****  |o(  >  +  X  (p.z)  \f,  > 

jl=-a>  i--® 


(B-1) 


where  (a)  and  (p)  are  the  "up"  and  "down"  spinors  respectively.  Apply  now 
and  Yq  to  Equation  (B-1).  Equations  (A-1),  (A-2)  and  Table  A-1  give: 


|1>  =  =  I  (p.z) 

£ 


ii((})±2ii/3)  „±1 


|a  > 


V  /  N  i£((bt2ii/3)  _  +  l  1-^ 

-h  X  (p.z)  e  '  n  |p> 


(B-2) 


B-1  Agassi,  D.  and  Restorff,  J.  B. ,  Pseudopotential  Band  Calculations  Along 
a  High-Svmmetrv  Axis:  Part  1- -Central  Potential  and  the  fill] 
Direction.  NAVSWC  TR  91-324,  15  Jun  1991. 
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|2>  -  n''^|2>  -  I  (p.z)  0"^  |a  > 

i 


V  17  /  \  ii(4)±2ii/3)  -+1 

+  X  (p,z)  e  '  '  Vi  |p> 


(B-3) 


Yq  |1>  -  |2>  -  X  (p.z)  e  |p>  -  X  (p.z)  e  |a  > 
i  i 


(B-4) 


where 


n  -  e 


(B-5) 


We  now  equate  the  spin  and  e^'^'^’  components  on  both  sides  of  Equations 
(B-2  -B-4).  Consider  first  Equation  (B-2).  Equating  the  U^(p,z)  components 
yields 


U^:  ^^1-^211/3  _  ^  =,  i  _  3m,  m  -  0,  ±1,  ±2, 


(B-6a) 


and  for  the  V^(p,z)  components 


V^;  0*^  =  1  =»  i  -  3m+l,  m  -  0,  ±1,  ±2, 


(B-6b) 


Similarly, 


(B-6d) 


W^:  0*^-1  =>  i  -  3m-l,  m  -  0,  ±1,  ±2,  ...  (B-6c) 

X^;  =  1  =>  i  -  3m,  m  -  0,  ±1 ,  ±2 ,  .  .  .  (B-Bd] 

To  represent  Equation  (B-6)  concisely,  we  introduce  the  following 
representation,  also  employed  in  Tables  3-1  and  3-2.  Since  |l>,  | 2>  in 
Equation  (B-1)  are  spinors 

^|a  >  -  ^  Q  j ,  |p  >  -  ^  ®  j  j,  the  U^(p,z),  V^(p,z)  and  W^(p,z),  X^(p,z) 

are  the  multipole  expansion  of  the  "up"  components  and  the  "down"  component, 
respectively.  Therefore  we  write  Equations  (B-6a)  and  (B-6b)  symbolically  as 


(B-7) 


11=--  [xial+i]  ■ 

I 1>  -  ^  3ni^+l  j  ^  shorthand  notation.  Equations  (B-6c)  and  (B-6d)  can 


similarly  be  written  as 
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1 


m  =  0,  ±1,  ±2,  ... 


(B-8) 


An  important  property  of  the  multipole  expansion  is  the  i-parity,  i.e., 
the  phase  relation  between  the  Jt  and  -i  multipoles.  This  information  is 
derived  from  Equation  (B-4): 

U_^(p.z)  =  X^(p,z) 


V_^(p,z)  =  -W^(p,z) 


(B-9) 


To  eliminate  the  p-dependence  from  Equation  (B-9)  recall  that,  e.g., 

Ujg(p.z)  =  Y  U^(g.z)  J^(gp)  (B-10) 


and  J_^(x)  =  (-1)'^  J^(x).  Therefore,  inserting  Equation  (B-10)  into 
Equation  (B-9)  gives 

(-1)^  U  ,(g,2)  =  X,(g,z) 


(B-lla) 


-(-1)"  V_^(g,z)  W^(g.z) 


(B-llb) 


where,  from  Equation  (B-2),  i  -  3m  or  i  -  3m+l.  Equation  (B-11)  express  the 
f -parity  relations.  They  can  be  combined  into  a  self  evident  s3nnbolic  matrix 
notation,  see  Table  3-3. 

The  other  two  spin  representations  of  are  analyzed  in  the  same 

manner.  For  the  double  group  we  take  the  equations  generated  by  X^.]^,  X^ , 

3d 

and  Yq  and  employ  Equations  (A-2)  through  (A-6)  and  Table  3-2.  The  results 
are  given  in  Table  3-4.  Note  that  a  representation  is  characterized  by 

the  z-parity,  i.e.,  the  phase  relation  between  e.g.,  U^(g,z)  and  U^(g,-z). 

A 

A  corrolar  result  is  obtained  by  comparing  the  .8 -parity  for  the 
allowed  i  values  Tables  3-1  and  3-2.  This  comparison  readily  yields  the 
"compatibility  relations"  between  representations  pertaining  to  a  A-point  and 
the  L-point®‘^: 


B-2  Slater,  J.  C.,  Quantum  Theory  of  Molecules  and  Solids.  Vol .  2,  McGraw- 
Hill,  New  York,  NY,  1965,  Appendix  9. 
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APPENDIX  C 

DERIVATION  OF  EQUATION  (4-3) 


Straightforward  algebra  leads  (see  Reference  C-1,  Appendix  C)  to  the 
following  expression 


-  (3  . 


Y  [  f(G')  X  (k'z'-f-L')]  e^^^z+^ 

G'  ,L' 


(C-1) 


In  Equation  (C-1)  z'  is  the  unit  vector  in  the  z'-axis  direction,  o  is  the 
Pauli  matrices  vector,  Equation  (2-8),  and  we  limit  ourselves  to  Ic'  vectors 
parallel  to  z'.  An  important  feature  of  Equation  (C-1)  is  that, 

aside  from  the  o  factor,  it  is  a  sum  of  products  of  two  factors  pertaining  to 

f(r')  and  ^  constrained  summation  and  z' -dependent  exponential. 

Therefore,  according  to  the  multipole  decomposition  theorem, 

expressed  in  terms  of  products  of  multipoles  pertaining  to  the  two  factors. 

To  carry  out  this  exercise,  we  first  focus  on  the  vector  product  in 
Equation  (C-1)  and  introduce  the  following  constructs: 


C-1  Agassi,  D.  and  Restorff,  J.  B. ,  Pseudopotential  Band  Calculations  Along 
a  Hi eh -Symmetry  Axis:  Part  I- -Central  Potential  and  the  fllll 
Direction.  NAVSWC  TR  91-324,  15  Jun  1991). 
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»<<’>, (gp.z’)  -  I  »<»>(£')  1. 

X ,  X  r  X 

-4 

L' 

ILtI-Sf 


»<“^,(gp,z')  .  1'  I  Ly,  . 

L' 

iL^hgp 


,i  v’  .Ao),:*,^  ^  ^i(L2.+k')z'  -ii6TlT') 


I  (k'-HL  .) 

A ,  Z  r  z 


(C-2) 


L' 

iL^I-gp 


By  noting  that  +  iLy  )/gF  =  e  T  (Reference  C-1),  these  multipoles 

take  the  form; 


(o 

i,x' 


(gp.z') 


>(o) 


(o), 


^^_y/(gF-z') 


>(o) 


(o), 


[^-(gp.z')  -’i'i:((gp.z')] 


9  . (o)  ,  , . 

az'  ^ 


(C-3) 


Each  components  of  the  vector  Equation  (C-1),  has  the  structure  of 

Equation  (C-1)  in  (Reference  C-1)  provided  v(2')  is  replaced  by  either 

fj^»(2'),  fyi(2')  or  f 2 > (S' ) (Equation  (2-14)),  and  ^(£')  is  replaced  by  either 

^(o)(j^')  T  ^  ^^*^^(1!')  h  or  )(k'+L_) .  Therefore,  the  decomposition 

theorem  (Appendix  C,  Reference  C-1)  applies  provided  Vjjj(gp,z')  is  replaced 

by  fjjj(gQ,z')  (Equation  (3-4)  in  main  text)  and  ^ ^  replaced  by 

Equation  (C-3).  These  steps,  and  the  evaluation  of  the  vector  product 

with  o  in  Equation  (C-1)  and  projection  on  both  (x(o))  spinors  gives 
Equation  (4-3)  of  the  main  text. 


» 
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APPENDIX  D 

MODELS  FOR  THE  SPIN-ORBIT  VECTOR  FORM  FACTOR 


The  spin  orbit  form  factor  f  (G')  is  related  to  the  radial  form  factor 
through  Equation  (2-11) 


f(G') 


-  I 
0  4 

J 


iG' ‘r i 
e  J 


7^  I/d?'  f  <|?'  -  ?;|)  1 

G'  0  J  J 


(D-1) 


In  Equation  (D-1)  0  is  ^he  unit  cell  volume,  j  runs  over  all  atoms  in  the  unit 
cell  with  coordinates  rj  with  respect  to  a  chosen  origin,  and  fj(|r'|)  is  a 
phenomenological  form  factor  for  the  j-th  atom. 


The  integral  in  Equation  (D-1)  should  be  carried  out  in  the  primed 
coordinates.  This  is  difficult  since  the  integration  limits  there  are 
awkward.  The  problem  is  avoided  by  transforming  the  integration  to  the 

unprimed  coordinated  and  then  applying  the  Vq*  operator.  In  the  unprimed 
coordinates  the  integral  in  Equation  (D-1)  takes  the  form 


I  (G)  =  /  dr  f  (|?-?J)  (D-2) 

J  Q  J  J 

Therefore,  for  a  separable  radial  form  factor,  1^(2)  takes  a  separable  form 
and,  hence,  easily  calculable.  The  expressions  below  for  a  constant  fj(|r- 
rj|)  are  easily  extended  to  a  separable  gaussian. 

For  a  constant  radial  form  factor  fj(|r-rj|)  =  Aj  with  a  range  r,  the 
integral  Equation  (D-2)  takes  the  form 

I.(G)  -  A.  J.(G  )  J.(G  )  J.(G  )  (D-3) 

J  Jjxjyjz 

where  the  integrals  depend  on  whether  the  atom  j  is  at  the  unit  cell 

center  (Se)  or  the  eight  eighth-atom  at  the  corners  of  the  unit  cell  (Pb). 

The  former  takes  the  form: 


D-1  Landsberg,  P.  T. ,  Solid  State  Theory.  Landsberg,  P.  T. ,  Ed.  Wiley- 
Interscience ,  London  1969,  p.  77. 
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-iG  t 

‘  ~  ^  ®  ""Se)  <D-^) 

■’^Se 

where 


2  sin  (G^r) 


G 

Tf 


Z2(0,A)  -  r 


(D-5) 


Consider  now  the  integrals  Ij(G^)  pertaining  to  the  Pb  atom.  For  an 
atom  at  the  r^  =  (0,0,0)  corner 

r 


0 


(D-6) 


By  the  same  token,  the  contribution  of  the  r2  -  (a, 0,0)  corner  is 


I 


(2) 

Pb 


a 

/  dx  e 


(x-a) 


a-r 


(D-7) 


Using  similar  relations,  the  Ip|j((5)  contributions  from  all  eight  corners  are 
easily  evaluated. 


We  turn  now  to  Equation  (D-1)  to  sum  the  contribution  from  all  atoms  in 
the  unit  cell.  The  tricky  parts  are  the  eight  Pb  contributions  located  at 


+  ?g  -  ?2  +  ^7  -  ^3  +  ^6  ■  ^4  ^5  "  ^  a*  (0,0,1)  (D-8) 


I 


and  a*  -  aj  3  .  With  the  help  of  Equation  (D-8),  the  form  factor 
Equation  (D-1)  takes  the  form: 


D-2 
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where 


exp  [  -  iG'-(?^  -  A'/2)  V-,  (l^CG^)  I^CGy)  ] 

-H  exp  [  -  iG’-(?^  -  AV2)  V-,  (I^(G^)  I^(Gy)  I^(G^))  ] 

+  exp  [  -  iG'-(?'  -  AV2)  vg,  (I^(G^)  I*(Gy)  ] 

-H  exp  [  -  iG'.(r;  -  AV2)  (I*(G^)  I*(Gy)  ]  | 

^Se  -  '^Sa  -■‘“'•''■''2  Sj.djCV  l2<V 


To  evaluate  Equation  (D-10),  it  is  necessary  to  carry  out  the  gradients. 
Since  the  gradient  transforms  like  a  vector  under  a  rotation,  we  can 
immediately  write: 


d 

A 

D 

d 

dc 

*  K 

ac 

y  J 

L  y 

(D-11) 


where  R  is  the  transformation  from  the  primed  to  the  unprimed  coordinates 
(Reference  D-2,  Equation  (2-1)).  Since  both  I2^(G.^)  and  12(G  )  are  given  in 
terms  of  elementary  functions,  Equation  (D-10)  can  be  evaluated  analytically. 


For  a  gaussian  form  factor,  the  above  analysis  remains  unchanged  except 
for  the  expressions  for  I2^(G^),  which,  instead  of  Equation  (D-5), 

involves  the  error  function. 


D-2  Agassi,  D.  and  Restorff,  J.  B.,  Pseudopotential  Band  Calculations  Along 
a  High-Svmmetrv  Axis:  Part  I --Central  Potential  and  the  fill] 
Direction.  NAVSWC  TR  91-324,  15  Jun  1991  (unpublished) 
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APPENDIX  E 

NUMERICAL  SOLUTION  OF  THE  MULTIPOLE  WAVE  EQUATION 


The  wave  equation  (Equation  (4-3)  of  the  main  text)  is  numerically 
solved  by  converting  it  to  a  secular  equation  as  described  in  Reference  E-1, 
Appendix  E.  We  elaborate  here  only  new  details  associated  with  the  SO 
interaction. 

To  ensure  that  the  numerically  computed  spin-orbit  form  factor 
(Equation  (2-11)  of  the  main  text)  has  the  correct  transformation  properties, 
we  invoke  the  analogy 


v'V(?')  <=>  f(r')  (E-1) 

which  follows  from  comparing  with  Equation  (2-8)  of  the  main  text.  In 
Equation  (E-1),  V(r')  transforms  as  a  scalar  under  the  groups  or  C^.^. 
Therefore,  from  Equation  (2-8) 


?  (G')  <->  i  G'  V(G') 


(E-2) 


or 


f  (G')  -  iG' 
av 


G'*f  (G') 
_ c _ 

[g'*g') 


(E-3) 


av 


In  Equation  (E-3),  f  (2')  denotes  the  calculated  spin-orbit  form  factor 


and 


G'-f  (G') 
_ c _ 

[g'*g'] 


Equation  (D-1)  and  is  the  average  over  all  2'  vectors 


av 


E-1  Agassi,  D.  and  Restorff,  J.  B. ,  Pseudopotential  Band  Calculations  Along 
a  High-Svmmetrv  Axis:  Part  I- -Central  Potential  and  the  fill] 

Direction .  NAVSWC  TR  91-324,  15  Jun  1991  (unpublished). 
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which  have  the  same  |2j|,  l^^l-  has  the  exact  transformation 

properties,  then  Otherwise,  has  the  correct 

transformation  properties.  The  symraetrization  Equation  (E-3)  has  been  used  in 
all  calculation. 


^m^  ^(gp,z')  or  to  the  derivative  d/dz' .  Formally,  the  derivative  appears  in 
a  hermitian  combination  (in  a);  f^"‘'^(z')  d/dz'  vs  -f^'^(z')  d/dz'  for  o-+l 

and  o=-l  equations,  respectively  (the  third  term  in  Equation  (4-3)).  This 
hermiticity,  however,  is  broken  in  the  process  of  converting  Equation  (4-3) 
into  a  secular  equation:  As  a  consequence  of  our  usage  of  sparse  Fourier 
components,  (which  depend  on  gQ  and  gp) ,  the  set  of  bra-plane  waves  |m>  may  be 
different  from  the  ket-plane  waves  <m]  (see  Appendix  E,  Reference  E-l)  whereas 
the  derivative  acts  always  to  the  right.  To  correct  for  this  artifact,  we 
symmetrize  the  derivative  term  as  follows  (Q(z)  is  an  arbitrary  function): 


The  o-mixing  terms  in  Equation  (4-3)  are  either  proportional  to 


(E-3) 


The  substitution  Equation  (E-3)  has  its  origin  in  the  spin-orbit  form, 

Equation  (2-8).  There  the  gradient  operator  acts  to  the  right.  Equivalently, 
the  gradient  can  be  made  act  to  the  left  by  commutation  with  the  spin-orbit 
form  factor  ,  i . e . , 


V  v(r)  X  (v)j^  -  -(v)j^  X  V  v(r)  (E-4) 

where  the  left  (L)  and  right  (R)  labels  explicitly  indicate  the  direction  of 
the  gradient  operator.  All  subsequent  calculations  employ  the  symmetrized 
derivative  form  Equation  (E-3). 
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